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Abstract 

For 117 clusters of galaxies, we explore the detectability of intracluster dust grains by current and 
future infrared facilities, taking into account both collisional heating and sputtering of grains by ambient 
plasma. If the dust grains are injected into the intergalactic space with the amount and size comparable 
to the Galactic values, the dust-to-gas ratio is typically 10 -6 and the mean dust temperature is ~ 30 K 
near the cluster center. The predicted infrared intensities lie marginally above the detection thresholds 
for Spitzer Space Telescope, ASTRO-F, Hcrschcl and SPICA missions. For some nearby clusters such as 
Perseus, A3571, A2319, A3112 and A2204, good detections of intracluster dust signal are expected in the 
70 fxxa band. Given rather tight constraints on the dust temperature from observed electron density and 
temperature, the dust mass can be inferred directly from the infrared observations. Further constraints 
on the size distribution will be obtained once multi-band data are available. They will definitely provide 
a powerful probe of the dust injection processes and dust-gas interactions in the intergalactic space. 

Key words: ISM: dust, extinction — galaxies: clusters: general — galaxies: intergalactic medium — 
infrared: general — X-rays: galaxies: clusters 



1. Introduction 

Dust grains have so far been detected only in the in- 
terstellar medium and there is no firm detection in the 
intergalactic medium. The presence of dust grains in 
the latter region is still an open question. If they exist, 
the intergalactic dust grains should have great impacts 
not only on our understanding of galaxy evolution but 
also on our interpretation of high redshift observational 
data. From X-ray measurements of heavy metal lines, it 
is evident that there is a significant amount of metal out- 
side of galaxies in clusters and groups of galaxies (e.g., 
Mushotzky et al. 1996; Rcnzini 1997; Davis, Mulchaey & 
Mushotzky 1999; Buote 2000). Such elements are likely to 
have been expelled out of galaxies through supernova ex- 
plosion, ram pressure stripping, or tidal interaction. Dust 
could also be expelled by these processes or by some inde- 
pendent mechanism such as radiation pressure (e.g., Chiao 
& Wickramasinghe 1972; Ferrara et al. 1991; Shustov & 
Vibe 1995; Davies et al. 1998; Simonsen & Hanncstad 
1999; Aguirre et al. 2001a,b) 

When dust grains arc placed in the intracluster medium 
(hereafter ICM), they would be heated by collisions with 
ambient hot X-ray emitting electrons and ions, and emit 
mainly in infrared bands (Dwek et al. 1990). At the 
same time, they are destructed via collisions with imping- 
ing ions (sputtering) with a typical timescale of ~ 10 8 
yr. The intracluster dust (hereafter ICD), if any, will 
therefore provide a direct clue to the ejection processes 
of materials from galaxies. Indeed, diffuse dust before de- 
struction has been detected in a number of elliptical galax- 
ies (e.g. Goodfrooij & de Jong 1995). The temperature, 



the amount, and the size distribution of the ICD, how- 
ever, can be quite different from those in the interstellar 
medium; in the absence of stellar photons, temperature 
of the ICD will be lower than that in elliptical galaxies. 
Moreover the ICD may be an important cooling agent in 
the ICM (Montier & Giard 2004). Revealing the nature of 
the ICD will be of particular importance in understanding 
the dust-gas interaction. 

There have been several suggestions and debates regard- 
ing the presence of the intergalactic dust. For example, 
Girardi et al. (1992) suggested that the redshift asym- 
metries of member galaxies in nearby groups are consis- 
tent with the presence of dust in the intragroup medium. 
The observed oxygen K edge in an X-ray spectrum of the 
Perseus cluster may also be attributed to the ICD (Arnaud 
& Mushotzky 1998). The extended submillimeter emis- 
sion detected in a rich galaxy cluster may partly be due 
to dust (Komatsu et al. 1999; Kitayama et al. 2004). On 
the other hand, the presence of enhanced visual extinc- 
tion towards high redshift objects behind clusters is still 
controversial (e.g., Maoz 1995). 

More direct evidence has been searched for in the far- 
infrared bands. Hickson et al. (1989) report that the 
far-infrared emission is enhanced by a factor of 2 in com- 
pact groups of galaxies compared with a sample of iso- 
lated galaxies. Sulentic & De Mello Rabaca (1993), how- 
ever, point out that the results of Hickson et al. (1989) 
arc likely to have been overestimated mainly due to the 
limited spatial resolution of IRAS. The emission statis- 
tically detected in the IRAS data with many clusters of 
galaxies may partly be due to the ICD (Montier & Giard 
2005). Stickel et al. (1998 and 2002) observed 6 galaxy 
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clusters with ISOPHOT, and only in Coma they reported 
the excess of ~ 0.2 MJy/sr at 120 /mi towards the cen- 
tral region. They estimated the dust color temperature of 
Id ~ 30 K, the excess flux of Ficd ~ 2.8 Jy and dust mass 
of M(j ~ 10 M @ from the ratio of surface brightness at 
120 /Ltm and 180 fxxa. The visual extinction derived from 
the dust mass is a negligible amount (Ay <C 0.1 mag). 
The result is in contrast to the reported optical extinc- 
tion of ~ 0.3 mag for Coma (Zwicky 1962, Karachentsev 
k Lipovetskii 1969). 

The major difficulties in the past infrared observations 
are the limited sensitivity and spatial resolution of the 
detector. The observational data thus might be affected 
by Galactic cirrus emission, contamination of individual 
galaxies, and confusion of extragalactic sources. The in- 
frared instruments on the current and future missions, 
such as Spitzer Space Telescope, ASTRO-F, Hcrschel, and 
SPICA, possess much improved sensitivity and spatial res- 
olution that are essential for separating the ICD emission 
from the other components. 

In this paper, wc explore the possibility of detecting the 
ICD with Spitzer Space Telescope, ASTRO-F, Herschel 
and SPICA by performing a comprehensive study of the 
expected intensities from the ICD for 117 clusters of galax- 
ies at redshift z ~ 0.01 — 0.8. Based on currently avail- 
able X-ray data, we explicitly take into account both 
heating and sputtering of dust grains via collisions with 
ambient hot plasma to compute their temperature and 
size distributions. We further examine the expected con- 
straints from the near future observations on underlying 
dust model. 

This paper is organized as follows. We describe our 
sample of galaxy clusters in Section 2. Details of the 
model for the ICD size distribution is presented in Section 
3. Section 4 describes the results for the properties of 
the ICD and the expected infrared spectra, as well as 
constraints from future observations on underlying dust 
model. We discuss several uncertainties in the present 
analysis in Section 5. Finally, Section 6 summarizes our 
conclusions. Throughout the paper, we assume a stan- 
dard set of cosmological parameters; f2 m = 0.3, ft a = 0.7, 
fi b = 0.04, and h = 0.7. 

2. Sample clusters 

In order to calculate the ICD emission, we pick up 117 
clusters at redshift z^O.01 — 0.8 for which gas properties 
can be determined by X-ray spectroscopy and imaging 
data from ROSAT and ASCA. Unless stated otherwise, 
we adopt an isothermal /3-model for the gas density profile 
of clusters. 

We extract 79 clusters at z ~ 0.1 — 0.8 from Ota & 
Mitsuda (2004). We also supplement the sample with 38 
nearby clusters at z ~ 0.01 — 0.1 from Mohr et al. (1999) 
and White (2000) for which the /3-model parameters, tem- 
perature, and metalicity are all available. For some clus- 
ters, Mohr et al. (1999) fit the X-ray surface brightness 
by a sum of two /3-model profiles. In such cases, we use 
the Abel's integral to deproject the fitted surface bright- 



ness profile into the electron density profile (Yoshikawa 
& Suto 1999). For the three clusters (A1689, A2204 and 
PKS0745-19) compiled in both Ota & Mitsuda (2004) and 
Mohr et al. (1999), we use the parameters from the former 
paper for definiteness. 

3. Model 

3.1. Dust- gas interaction 

We first need to specify the basic dust model to calcu- 
late the ICD emission. We consider a mixture of equal 
amounts of spherical silicate and graphite as dust com- 
ponents, and adopt as their mass densities Pgra — 

2.26 g 

cm" 3 (Draine & Lee 1984) and p s \\ = 3.5 g cm" 3 (Li & 
Draine 2001), respectively. Using these two components, 
the extinction curve in our Galaxy can be well repro- 
duced (Mathis, Rumpl & Nordsieck 1977, hereafter MRN; 
Draine & Lee 1984; Laor & Draine 1993). 

If dust grains exist in the ICM, they will be heated by 
collision with ambient X-ray emitting hot electrons and 
ions, and consequently emit in infrared band as they cool. 
The temperature of the ICD is then obtained via the bal- 
ance between collisional heating and radiative cooling. 

For given gas temperature and density, we determine 
the temperatures of dust grains following Dwek (1986) 
taking account of both heating and cooling between each 
collision. The heat capacities for graphite and silicate are 
taken from Dwek (1986) and Draine & Anderson (1985), 
respectively. Figure 1 shows a specific example for the 
temperature distribution of graphite grains embedded in a 
plasma with electron temperature T c = 10 s K and electron 
density n c = 10 -3 cm" 3 . Note that the dust temperature 
is not uniquely specified for given electron temperature 
and density, but has a probability distribution depending 
on the grain size. For small (a < 0.1/im) grains, the cool- 
ing time gets shorter than the interval of collisions and the 
dust temperature has a large dispersion (so-called stochas- 
tic heating). Some portion of these grains can have rather 
high temperatures ~ 100 K owing to their small heat ca- 
pacity. Larger grains, on the other hand, approach ther- 
mal equilibrium with ambient plasma and the tempera- 
ture distribution has a sharp peak at ~ 20 K. As a result 
of stochastic heating, the spectral shape of dust emission 
will be broadened particularly in the Wien regime. We 
will explicitly take into account such features in our pre- 
diction of the infrared emission from the ICD. 

When dust grains are placed in a hot ICM, they are 
sputtered by impinging ions. We adopt the following ana- 
lytic form for the sputtering rate (Tsai & Mathews 1995): 



da 
~~dt 



-r- = ~h 



Pgas 



T 



1 



(1) 



which gives a good approximation to the detailed calcu- 
lations of Draine & Salpeter (1979) for both graphite and 
silicate grains when ft = 3.2x 10 -18 cm 4 s _1 , T s = 2x 10 6 
K, and w = 2.5. Using this formula, wc define the sput- 
tering timescale as 



T sp utt(a,r) = a \ da/dt 



(2) 
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Fig. 1. Temperature distribution of graphite grains with var- 
ious sizes embedded in a hot plasma with electron density 
10 — 3 cm -3 and electron temperature 10 s K. 

For gas temperatures above ~3x 10 6 K, the sputtering 
timescale is almost independent of gas temperature and 
approximated by 



T spu tt(a,r) ~ 10 8 



io- 



yr. 



(3) 



0.1/xm 

3.2. Size distribution of dust grains in the ICM 

In modelling the size distribution of the ICD, we assume 
that the dust grains are continuously injected into the 
ICM from galaxies and destructed via sputtering. The 
size distribution dni(a;r)/da for species i at radial point 
r then follows (e.g., Dwek et al. 1990; Laor & Draine 
1993) 



d_ 
da 



. 9rii(a;r)^ | d dni(a;r) 9n™ J (a;r) 



do 



dt da 



da 



cluster 



.,(4) 



where d = da/dt is given by equation (1), dn i i {a;r)/da 
describes the total number density of injected dust grains 
in the size interval between a and a + Sa at radial point r, 
and T c i us ter is the cluster lifetime over which dust grains 
have been injected. Neglecting uncertainties for the for- 
mation epochs of individual galaxies in clusters, we set 
^cluster equal to the age of the Universe at the cluster red- 
shift. For the injected dust grains, we assume that they 
trace the galaxy distribution in a cluster and have the size 
distribution observed in our Galaxy (MRN): 



(9n- nj (a;r) 
da 



= A l {r)a~ 



for a„ 



< a < a r 



(5) 



where a lnj = 3.5, a m - m = 0.001 /xm, a max = 0.25 /jm, and 
Ai(r) oc p ga ,i(r) is the normalization coefficient to be de- 
termined later. Equation (4) is solved for steady-state to 
give 



dnf d (a;r) _ 

d~a ~ Mr) 



a T sput t(a,r) 
a inj - 1 T cluster 



(6) 



Since r spu tt oc a, the steady-state size distribution is pro- 
portional to a _Q ' nJ+1 . 

Given large uncertainties in the galaxy distribution in 
each cluster, particularly at its outer envelope, we assume 
that galaxies trace collisionless dark matter and Ai(r) oc 
Pgal(^) oc /9 DM (r). We adopt the density profile of dark 
matter inferred from N-body simulations (Navarro, Frenk 
& White 1996, hereafter NFW): 

PdmM oc - 1 — (7) 

{r/r s ){l + r/r s y 

where r s is the scale radius of the halo. We follow Bullock 
et al. (2001) to determine r s for a halo with virial mass 
M v ; r at redshift z. As in Shimizu et al. (2003), we trans- 
late iteratively the scaled mass M500 , given for each cluster 
in Mohr et al. (1999) and Ota & Mitsuda (2004), into M vir 
using the NFW profile, where M500 is the mass enclosed 
within the radius at which the mean density p(r) is equal 
to 500 times the critical density of the universe at the 
cluster redshift. The dark matter density is normalized 
by M vir = f rvir p UM (r)4irr 2 dr. 

Taken together, Aj(r) is determined by fixing the total 
amount of the injected dust grains M^ nj as 



AC(<rvir) = 



E 



4irr 2 Ai(r)dr 



-,3-a" 



da 



— Z d J ^metalM gas (< r v j r ), 



(8) 



where the injected dust-to-gas mass ratio Z™ s is fiducially 
fixed at 0.0075 (MRN), and the mean metalicity Z meta i 
in units of the solar value is taken from Ota & Mitsuda 
(2004) and White et al. (2000). We have assumed that 
ZmetaiM gas (< r v ; r ) denotes the cumulative mass of the gas 
that has been supplied from galaxies into the ICM. Using 
the assumed relation of Ai (r) oc p DM (r) , we obtain 



Mr) 



3/i Z A Z meta \ -Mg as (< r v ; r ) 



PdmW, (9) 



An Pi f^a 3 - aini da 

where fa represents the mass fraction of graphite or silicate 
grains. 

3. 3. Infrared emission from the ICD 

Given the size distribution (eq. [6] ) and the temperature 
distribution Gi(a,T^) of the ICD, the infrared intensity at 
wavelength A is given by an integral over the line-of-sight 
through the cluster: 




da na Q\(a) 



Gi{a,T A )B x {T d )dT d 



dnf d (a;r) 



da 



(10) 



where Q\(a) is the dust absorption efficiency factor of 
dust species i taken from Draine & Lee (1984), B\ is the 
Planck function, R is the projected separation from the 
cluster c enter, and I is the length over a line of sight with 
'max = \j '"vir — R 2 • ^hc minimum temperature of the dust 
Tmi n is set equal to the CMB temperature 2.73(1 + z) K, 
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and the maximum temperature T max to the vaporization 
temperature, 2000 K for graphite and 1500 K for silicate, 
respectively (Dwek 1986). As Gi(a,T d ) declines rapidly 
toward high temperatures, our results are insensitive to a 
specific choice of T max . 

4. Results 

4-1. Properties of the ICD 

Based on the model constructed in the previous section, 
the ICD mass density as a function of physical radius r 
from the cluster center is given by 

lax Ana 3 pi dnf d (a;r) 



Pd 



M = £/< 



da 



-da. 



(11) 



The dust-to-gas mass ratio in the ICM can then be written 

as 

Pd(r) 



Mr) 



(r) 



(12) 



We define the mass weighted mean temperature and mean 
size of the ICD as follows: 



<Td(r)> = 




Gi(a,T d )T d dT d 



Ana 3 Pi dnf d (a:r) , 
x l — - -da 



Mr)) 



da 
4na 3 pi 



Pd(r), 



(13) 



dnf d (a:r) , 
x — 1 ^ ada 
da 



Pd(r) 



(14) 



For later convenience, we further introduce the surface 
mass density of the ICD at a projected separation R from 
the cluster center by 

Ed(#) = J ^ p d (r)dl. (15) 

Figures 2-4 show the above quantities in a specific case 
of Perseus cluster (A426). As will be shown below, the 
expected intensity from this cluster is among the highest 
in our entire sample and we will take it as a representative 
target in our subsequent analysis. The dust-to-gas mass 
ratio Z d (r) is lower, particularly near the center, by sev- 
eral orders of magnitude than the Galactic value 0.0075 as 
a result of rapid sputtering (Fig. 2). The mean dust tem- 
perature (Td(r)), on the other hand, is the highest near the 
center at > 30 K due to efficient heating by electron colli- 
sions (Fig. 3). The mean size of the dust is independent of 
r and its value in Perseus is ~ 0.15 /xm. Once integrated 
over the line-of-sight, the ICD surface mass density Ed 
remains nearly constant within the projected separation 
of ~ 300 kpc (Fig. 4) and will serve as a good measure of 
the total amount of the ICD. 

4-2. Expected infrared spectra 

We plot the predicted ICD emission spectra from 
Perseus cluster in Figure 5. In the rest of the paper, 
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Fig. 2. Dust-to-gas mass ratio in the ICM of Perseus cluster 
(A426) as a function of physical radius from the center. 




Fig. 3. Mass weighted mean temperature of the ICD in 
Perseus as a function of physical radius from the center. 




10 1 
R [Mpc] 



Fig. 4. Surface mass density of the ICD in Perseus as a func- 
tion of projected separation from the cluster center. 
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Fig. 5. Predicted spectra of the ICD at the projected sepa- 
ration of 20 kpc from the center of Perseus cluster. Open cir- 
cles, open triangles and open squares show the 3cr confusion 
limits, due to extragalactic sources, for ASTRO-F, Spitzcr, 
Herschel and SPICA missions, respectively. Filled symbols 
with an identical shape indicate those due to Galactic cir- 
rus for the same instrument. Note that the open and filled 
triangles overlap each other at 160 (im. 

unless otherwise stated explicitly, the emission is com- 
puted at the projected separation of 20 kpc from the cen- 
ter in order to avoid any contamination from the central 
galaxy. Also plotted for reference are the 3cr confusion lim- 
its due to extragalactic sources, for Spitzer/MIPS (Dole et 
al. 2004), ASTRO-F/FIS (Pearson ct al. 2004), Herschel 
and SPICA (Dole et al. 2004), as well as those due to 
Galactic cirrus (Jeong et al. 2005). The confusion levels 
are computed at the beam size of each instrument. For 
simplicity, we represent band filters of FIS (N60, N170, 
WIDE-S and WIDE-L) by their central wavelengths 60, 
170, 75 and 150 /im, respectively. We adopt 70 /im and 
160 (im for the band filters of Herschel and SPICA. The 
PSF of each instrument is modeled by a Gaussian with 
the FWHM value taken from the reference mentioned 
above. We also take into account the dependence of the 
cirrus emission on the Galactic latitude as in Jeong et al. 
(2005). The total confusion limit should be estimated as 

CT tot = ^extragalactic + (J cirru S - This n S Ure indicates that 

the source confusion is dominated by extragalactic sources 
at A < 100/iin, while the contribution of Galactic cirrus 
increases rapidly at longer wavelengths. The predicted 
emission has its peak at ^ 70— 100 /im and would be de- 
tectable above 5cr levels at 70 /im with Spitzer, Herschel 
and SPICA, and marginally detectable at the 3c levels at 
60 and 75 /im with ASTRO-F. 

Figure 6 shows the expected radial profiles for the emis- 
sion from Perseus. The emission at 70 /zm would be de- 
tectable above the 3er levels within ~ 2.5 arcmin by Spitzcr 
and within ~ 4.5 arcmin by Herschel and SPICA. At the 
outer envelopes of the cluster, the emission declines owing 
to the lack of hot gas capable of heating the dust. The 
extended feature of the emission will be of particular im- 
portance in separating the ICD component from galaxies. 



Fig. 6. Expected radial profiles of the ICD emission at 70 fira 
(thick solid line) and 160 (im (thin solid line) in Perseus clus- 
ter. Also plotted for reference are the 3<r confusion limits at 70 
£im; extragalactic sources for Spitzer (dashed line), Herschel 
and SPICA (dash-dotted line), and Galactic cirrus for Spitzcr 
(dotted line), Herschel and SPICA (dash-dot-dotted line). 

A fine spatial resolution of Spitzcr, ASTRO-F, Herschel 
and SPICA will be essential for this purpose. 

In Figure 7, we plot the expected intensities of the 
ICD for the entire sample of clusters as a function of the 
Galactic latitude. In general, the highest signal-to-noise 
ratio (S/N) is achieved at 70 /im. Five clusters with the 
highest 70 /im intensity (filled symbols) all lie above the 4cr 
levels for Spitzer and above the 7er levels for Herschel and 
SPICA. The emission tends to decrease toward high red- 
shift because the observed intensity and wavelength vary 
as oc (1 + z)~ A and oc 1 + z, respectively. Several nearby 
clusters at z < 0.1 lie above the 3a levels of the detector at 
70 /zm, 150 /tm, and 170 /im. They are likely to serve as 
plausible candidates for the positive detection of the ICD. 

We note that the confusion noise per beam, due to ex- 
tragalactic sources, drops slower than the size of the aper- 
ture squared at A > 100 /im (e.g., Dole et al. 2004). This 
makes the extragalactic confusion per fixed sky area in- 
crease with an increasing telescope size (Figs 5 and 7 ). 
At such long wavelengths, one may be able to achieve 
higher S/N by adding together the larger sky area, while 
high resolution detectors should still be useful to exclude 
individual galaxies. 

In order to examine the scale dependence of the ICD de- 
tectability, we plot in Figure 8 the average surface bright- 
ness of A3571 at 160 /im within a given radius excluding 
the central 20 kpc region. Also plotted for reference are 
the confusion levels for the same area. We have chosen 
A3571 because it has the highest S/N at 160 /zm among 
the five representative clusters mentioned above. In com- 
puting the confusion noises for larger sky area, we have 
extrapolated the values indicated in Table 1 and 2 of Dole 
et al. (2004) for extragalactic sources and Figure 17 of 
Jeong et al. (2005) for Galactic cirrus, respectively. The 
extragalactic noise drops more rapidly with scale than 
the cluster signal, while the Galactic confusion increases 
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Fig. 7. Expected intensities at 60, 70, 75, 150, 160 and 170 (im at the projected separation of 20 kpc from the center of 117 galaxy 
clusters versus Galactic latitude. Five clusters with the highest 70 (im intensities are marked by filled symbols; Perseus (circle), 
A2319 (square), A3571 (lozenge), A2204 (triangle), and A3112 (reverse-triangle). The other clusters are marked by open circles; 
at redshift 0.01 < z < 0.1 (large circles), at 0.1 < z < 0.3 (medium circles), and at 0.3 < z < 0.8(small circles). Lines represent 3<r 
confusion limits, due to extragalactic sources and Galactic cirrus for the mission indicated in each panel. 
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Fig. 8. Average surface brightness of A3571 at 160 
within a given radius excluding the central 20 kpc region. 
Also plotted for reference are the confusion levels at 160 fim 
for the same area; due to extragalactic sources (dashed line) 
and Galactic cirrus (dotted line). 

monotonically. The highest S/N is expected within 6* ~ 45 
arcsec. This sort of multi-scale analysis is a powerful tool 
for detecting the ICD at > 100 /zm. 

4-3. Constraints on the ICD properties 

We have so far assumed that the dust has been injected 
into the ICM with essentially the same size distribution 
and amount as in our Galaxy. We now examine how sensi- 
tive our results are to their uncertainties and explore what 
constraints one can place in turn on the underlying dust 
model from future observations. 

Specifically, we vary three parameters in our model; the 
power-law index a of the dust size distribution, their max- 
imum size a max , and the surface mass density Ed at the 
projected radius R = 20 kpc. Since very small dust grains 
are rapidly destroyed in the ICM and do not contribute to 
mid- to far-infrared emissions, we simply fix the minimum 
dust size at a m i n = 0.001 /zm. As mentioned in Section 3, 
the power-law index a is related to that of the injected 
dust with a = a" 11 — 1 and has the value 2.5 for a lns = 3.5 
(MRN). As shown in Figure 4, Ed is nearly constant in the 
cluster and serves as a good measure for the total amount 
of the ICD. We take Perseus as a representative target and 
use the data at 24, 70 and 160 /zm to derive constraints 
on the three parameters mentioned above. 

Figure 9 shows contours of the 70 /urn intensity on 
«max — Ed plane for a = 2.5 and a — Ed plane for et max = 
0.25 /zm. We find that the intensity is much more sensitive 
to Ed than to a max and a. This indicates that the detec- 
tion of (or an upper limit to) the ICD emission directly 
leads to a measure of the total amount of dust almost inde- 
pendently of its size distribution. Note that the mean ICD 
temperature varies only moderately and can be well con- 
strained once electron density and temperature are known 
from X-ray observations (Fig. 3). 

Constraints on the size distribution, a max and a, will 
be obtained via multi-band observations as illustrated in 



Figure 10. Since the intensities are proportional to Ed, 
uncertainties in Ed are eliminated by taking their ratios 
in multi-bands. The ratio 7(160 /zm)/7(70 /zm) can yield 
a reasonable measure of a max provided that a is not sig- 
nificantly different from that inferred in the MRN model 
a = 2.5. On the other hand, 7(24/zm)/7(70 /zm) is more 
sensitive to a than a max and can be used to infer the for- 
mer quantity if sufficient S/N is achieved in both bands. 

5. Discussion 

There have been a number of suggestions as to how 
dust grains arc replenished into the intergalactic medium. 
For example, many authors have argued that a signifi- 
cant amount of dust could be expelled from galaxies by 
radiation pressure (e.g., Chiao & Wickramasinghe 1972; 
Ferrara et al. 1991; Shustov & Vibe 1995; Davies et 
al. 1998; Simonsen & Hannestad 1999; Aguirre et al. 
2001a, b). Radiation pressure may limit the injected dust 
size at 0.05 — 0.2 /zm (Shustov & Vibe 1995 ; Davies et 
al. 1998) because smaller grains are rapidly destroyed 
via sputtering in galactic halos while larger ones are hard 
to escape from the gravitational potential of galaxies. 
Ejection of small grains can be suppressed further by drag 
forces due to dust-gas collisions and Coulomb interactions 
(Bianchi & Ferrara 2005). Ferrara et al. (1991) point out 
that radiation pressure may also change the composition 
of ejected grains because graphite grains can attain higher 
velocities and escape more easily from galaxies than sili- 
cate grains. 

Another potential source of the ICD is from the intra- 
clustcr stellar population (Montier & Giard 2004). Zwicky 
(1951) was the first to suggest the presence of the intra- 
cluster stars based on the detection of the excess light 
between the galaxies of the Coma cluster. Numerical sim- 
ulations indicate that ~ 10 — 20 per cent of the total cluster 
light should be the ICM (Murante et al. 2004; Willmanet 
et al. 2004). Recently, planetary nebula (e.g. Feldmeier et 
al. 2004; Theuns & Warren 1997), red-giant-branch stars 
(Ferguson et al. 1998; Durrcll et al. 2002), and also Type 
la supernova (Gal- Yam et al. 2003) have been detected in 
the ICM. The amount of intracluster light allowed from 
these observations is ~ 10 — 50 per cent of the total cluster 
light. According to Montier & Giard (2004), the dust-to- 
gas mass ratio of the Coma-like cluster, due to the intra- 
cluster stellar population, is ~ 10~ 5 within a radius of 1 
Mpc and predominant compared to the galactic produc- 
tion. The intracluster stellar population can add small 
dust and silicate grains that might be suppressed if galax- 
ies are the only source of the ICD. As demonstrated in 
Section 4.3, multi-band observations of the dust emission 
will provide a powerful probe of the grain size distribution, 
which can be used to constrain the ejection processes of 
the ICD. 

6. Conclusions 

In this paper, we have performed a comprehensive study 
on the nature of diffuse dust in a sample of galaxy clusters 
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Fig. 9. Contours of the 70 fim intensity from Perseus on an a max — plane for a = 2.5 (left panel) and an a — plane for 
fimax = 0.25/im (right panel). The curves correspond to 7(70 (im) = 100, 10 and 1 mJy/arcmin 2 from top to bottom in the both 
panels. 




Fig. 10. Contours of the intensity ratios, 7(24 /xm)/7(70 fim) (left panel) and 7(160 fim)/ 1(70 fim) (right panel) on an a max — ot 
plane for Perseus. The curves, from top to bottom, correspond to 7(24 fim)/I(70 fim) = 0.1, 0.01 and 0.001 in the left panel and to 
7(160 £*m)/7(70 ^m) = 0.4, 0.6, 0.8 and 1.0 in the right panel. 
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at redshift z ~ 0.01 — 0.8. Based upon recent X-ray data, 
the temperature distribution of the dust grains are com- 
puted taking account of collisional heating by ambient hot 
plasma and radiative cooling. The dust size distribution 
is also solved by incorporating injection from galaxies and 
destruction via sputtering. 

If the dust grains are injected into the intergalactic 
medium with the amount and size comparable to the 
Galactic values, the dust-to-gas ratio is typically ~ 10~ 6 
and the mean dust temperature is ~ 30 K near the cluster 
center. The predicted emissions lie marginally above the 
detection thresholds for Spitzer Space Telescope, ASTRO- 
F, Herschel and SPICA missions. The highest signal-to- 
noise ratios (> 4cr) are expected for some nearby clusters 
such as Perseus, A3571, A2319, A3112 and A2204 at the 
70 band. 

Given rather tight constraints on the dust temperature 
from observed electron density and temperature, the dust 
mass can be inferred directly from the infrared observa- 
tions. Further constraints on the size distribution will be 
obtained once multi-band data are available by the future 
facilities. They will definitely provide a powerful probe of 
the dust injection processes and dust-gas interactions in 
the intergalactic space. 

We thank Woong-Seob Jeong for offering the confusion 
noise data due to Galactic cirrus, Naomi Ota for providing 
the X-ray data of distant clusters, and the referee Martin 
Giard for helpful comments. We also thank Hidehiro 
Kaneda, Hirohisa Nagata, Mamoru Shimizu, and Hidenori 
Takahashi for useful discussions. This work is supported 
in part by the Grants- in- Aid by the Ministry of Education, 
Science and Culture of Japan (14740133:TK). 
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